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The Dynamical Cluster Approximation (DCA) is modi- 
fied to include disorder. The DCA incorporates non-local 
corrections to local approximations such as the Coherent 
Potential Approximation (CPA) by mapping the lattice 
problem with disorder, and in the thermodynamic limit, 
to a self-consistently embedded finite-sized cluster problem. 
It satisfies all of the characteristics of a successful cluster 
approximation. It is causal, preserves the point-group and 
translational symmetry of the original lattice, recovers the 
CPA when the cluster size equals one, and becomes exact as 
7V C — *• oo. We use the DCA to study the Anderson model 
with binary diagonal disorder. It restores sharp features 
and band tailing in the density of states which reflect cor- 
relations in the local environment of each site. While the 
DCA does not describe the localization transition, it does 
describe precursor effects of localization. 



I. INTRODUCTION 

The Coherent Potential Approximation (CPA) 
is a widely used method for treating disordered sys- 
tems. Within the CPA, the problem is first averaged 
over all possible disorder configurations in an attempt 
to regain the translational invariance lost due to dis- 
order. Then non-local correlations of the disorder po- 
tential are neglected leading to a self consistent single- 
site approximation (like the Weiss mean field theory of 
magnetism). It has been applied with great success to 
a variety of problems , including ab-initio calcula- 
tions in disordered metallic alloys Q]. 

Nevertheless, the CPA fails to provide a completely 
satisfactory theory for disordered systems ||. As a 
single-site mean field theory, it cannot account for the 
disorder induced, short ranged but nonlocal, correla- 
tions due to the local order in the environment of each 
site responsible for band tailing and sharp structures in 
the density of states. There have been many attempts 
HHI to formulate non-local corrections to the CPA in 
which the lattice is mapped onto a self-consistently em- 
bedded finite-sized cluster. However, as argued clearly 
and in detail by Gonis j|, these theories all fail in 
some significant way. A successful theory must be able 
to account for fluctuations in the local environment 
in a self-consistent way, become exact in the limit of 
large cluster sizes, and recover the CPA when the clus- 
ter size equals one. It must be easily implementablc 



numerically and preserve the translational and point- 
group symmetries of the lattice. Finally, and most sig- 
nificantly, it should be fully causal so that the single- 
particle Green function and self energy are analytic in 
the upper half plane. No presently existing cluster ex- 
tension of the CPA satisfies all these requirements ||,|7j . 

Recently, a new method called the The Dynamical 
Cluster Approximation (DCA) |^— [To|] was developed 
for ordered correlated systems such as the Hubbard 
model to add non-local corrections to the Dynamical 
Mean Field Approximation. In this manuscript, we 
modify the DCA to include disorder and show that 
the resulting formalism satisfies all of the above stated 
requirements for a successful cluster extension of the 
CPA. 

In the next Section, we review the basic diagram- 
matic perturbat ion t heory formalism for disordered 
systems. In Sec. Ill A we show that the CPA is equiva- 
lent to neglecting momentum conservati on at all inter- 
nal vertices of the diagrams, and in Sec. IIIB we intro- 
duce the DCA for disordered systems which systemati- 
cally restores the momentum conservation relinquished 
by the CPA. In Sec. |^ we show that the DCA satis- 
fies each of the desired characteristics described above. 
In Sec. we show results for the two-dimensional An- 
derson model with binary diagonal disorder. Finally in 
Appendix A we present an alternate way of viewing and 
justifying the disordered DCA algorithm developed in 
this paper, using the replica (or other) methods of dis- 
order averaging. 



II. BASIC FORMALISM 

We consider an Anderson model with diagonal dis- 
order, described by the Hamiltonian 

H = - E 1 ( C L C ^ + C L C ^) + - a*k„ 



<y>,c 



(1) 



where C\ a creates a quasiparticle on site i with spin 

= C\ a C iitr . The disorder occurs in the lo- 
cal orbital energies Vi, which we assume are indepen- 
dent quenched random variables distributed according 
to some specified probability distribution P(V). The 
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DCA formalism that we develop in this paper is a gen- 
eral method valid for any P(V). However, for illus- 
trative purposes, for the specific calculations presented 
in this paper we take V* = ±V with equal probability 
1/2 (binary disorder). 

The effect of the disorder potential Ylicr Vi n i,<? can be 
described using standard diagrammatic perturbation 
theory (although we will eventually sum to all orders), 
ft may be re-written in reciprocal space as 
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(2) 



The corresponding irreducible (skeletal) contributions 
to the self energy may be represented diagrammatically 
H and are displayed in Fig. |l|. 



+ 



5 2 



+ 



+ 



13 4 5 2 1 3 4 5 6 2 

FIG. 1. The first few graphs in the irreducible self en- 
ergy of a diagonally disordered system. Each o represents 
the scattering of a state k from sites (marked X) with a 
local disorder potential distributed according to some speci- 
fied probability distribution P(V). The numbers label the k 
states. 

Here each o represents the scattering of an electronic 
Bloch state from a local disorder potential at some site 
represented by X. The dashed lines connect scattering 
events that involve the same local potential. In each 
graph, the sums over the sites are restricted so that 
the different X 's represent scattering from different 
sites. No graphs representing a single scattering event 
are included since these may simply be absorbed as a 
renormalization of the chemical potential /i. 

Translational invariance and momentum conserva- 
tion are restored by averaging over all possible values 
of the disorder potentials Vi. For example, consider the 
second diagram in Fig. [I], given by 



^3 (^ 3 }G(k 3 )G(k 4 ) e ^( kl -k 3+k3 - 

i,k 3 ,k 4 



-k 4 +k 4 — k 2 ) 



(3) 



where G(k) is the disorder-averaged single-particle 
Green function for state k. The average over the distri- 
bution of scattering potentials (V?) = (V 3 ) indepen- 
dent of i. After summation over the remaining labels, 
this becomes 



(V 3 )G(r 



0) 2 S kl ,k 2 



(4) 



where G(r = 0) is the local Green function. Thus 
the second diagram's contribution to the self energy 
involves only local correlations. Since the internal mo- 
mentum labels always cancel in the exponential, the 
same is true for all non-crossing diagrams shown in the 
top half of Fig. 0. 

Only the diagrams with crossing dashed lines are 
non-local. Consider the fourth-order diagrams such 
as those shown on the bottom left and upper right of 
Fig. [|. When we impurity average, we generate po- 
tential terms (V 4 ) when the scattering occurs from the 
same local potential (i.e. the third diagram) or (V 2 ) 2 
when the scattering occurs from different sites, as in 
the fourth diagram. When the latter diagram is evalu- 
ated, to avoid overcounting, we need to subtract a term 
proportional to (V 2 ) 2 but corresponding to scattering 
from the same site. This term is needed to account 
for the fact that the fourth diagram should really only 
be evaluated for sites i ^ j\ For example, the fourth 
diagram yields 



_ \ " V 2 V 2 p iMki+k 4 -k 5 -k 3 ) ir,-(k 5 +k 3 -k 4 -k 2 ) 

G(k 5 )G(k 4 )G(k 3 )) (5) 
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Evaluating the disorder average (), we get the following 
two terms: 

_L ^2^2 e ir ! -(ki+k 4 -k5-k 3 ) e ir j -(k 5 +k 3 -k 4 -k 2 ) 
ijk 3 k 4 k 5 

G(k 5 )G(k 4 )G(k 3 ) 
-Jji E (^ 2 ) 2 e^ (kl - k2) G(k 5 )G(k 4 )G(k 3 ) (6) 

ik 3 k 4 ks 

Momentum conservation is restored by the sum over i 
and j; i.e. over all possible locations of the two scatter- 
ers. It is reflected by the Laue functions, A = iV<5k+-.-, 
within the sums 



^k 2 ,ki ^ 
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G(k 5 )G(k 4 )G(k 3 ) 
(F 2 ) 2 G(k 5 )G(k 4 )G(k 3 ) 



(7) 



Since the first term in Eq. |?] involves convolutions of 
G(k) it reflects non-local correlations. [Local contri- 
butions such as the second term in Eq. M can, if one 
so chooses, be combined together with the contribu- 
tions from the corresponding local diagrams such as 
the third diagram in Fig. p] by replacing {V 4 ) in the 
latter by the cumulant (V 4 ) — (V 2 ) 2 .} Given the fact 
that different X's must correspond to different sites, it 
is easy to see that all crossing diagrams must involve 
non-local correlations. 
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III. CLUSTER APPROXIMATIONS 



B. The Dynamical Cluster Approximation 



A. The Coherent Potential Approximation 

In the Coherent Potential Approximation (CPA), 
nonlocal correlations involving different scatterers are 
ignored. Thus, in the calculation of the self energy, we 
ignore all of the crossing diagrams shown on the bot- 
tom of Fig. [|; and retain only the class of diagrams 
such as those shown on the top representing scattering 
from a single local disorder potential. These diagrams 
are shown in Fig. 0. 




FIG. 2. The first few graphs of the CPA self energy. Here 
the Green function line represents the average local propa- 
gator. 



Employing diagrammatic arguments, it is easy to see 
that the CPA is fully equivalent to the neglect of mo- 
mentum conservation at each internal vertex. This is 
accomplished by setting each Laue function within the 
sum (eg., in Eq. [7]) to 1. We may then freely sum 
over the internal momenta, leaving only local propaga- 
tors. All non-local self energy contributions (crossing 
diagrams) must then vanish. For example, consider 
again the fourth graph. If we replace the Laue func- 
tion A r <5k 1 +k 4 .k5+k3 1 in Eq. [j], then the two con- 
tributions cancel and this diagram vanishes. Thus an 
alternate definition of the CPA (TTJ , in terms of the 
Laue functions A, is 

A = Acpa = 1 (8) 

I.e., the CPA is equivalent to the neglect of momen- 
tum conservation at all internal vertices of the disorder- 
averaged irreducible graphs. 




FIG. 3. N c — 4 cluster cells (shown by different fill pat- 
terns) that partition the first Brillouin Zone (dashed line). 
Each cell is centered on a cluster momentum K (filled cir- 
cles). To construct the DCA cluster, we map a generic 
momentum in the zone such as k to the nearest cluster 
point K = M(k) so that k = k — K remains in the cell 
around K. 

The DCA systematically restores the momentum 
conservation at internal vertices which was relin- 
quished by the CPA, and so incorporates non-local cor- 
rections. This is done by dividing the Brillouin zone 
into N c equal cells, as shown in Fig. || and requiring 
that momentum be partially conserved for momentum 
transfers between the coarse graining cells shown in 
Fig. ^, but ignored for momentum transfers within each 
cell. This may be accomplished by employing the Laue 
functions 

A = A DC A = -^c<5M(k 1 )+M(k 2 ),M(k 3 )+M(k 4 )... (9) 

where M(k), illustrated in Fig ||, maps the momenta 
k to the nearest cluster momenta K. Adca becomes 
one when iV c = I since then all momenta are mapped 
to the zone center by M. Thus the CPA is recovered 
in this limit. Furthermore, as iV c becomes large, the 
exact result is recovered since limjv c ^oo M(k) = k for 
all momenta k. 

>, A DCA = Nc§ M(kl) + M(k2) _ M(k3) + M(k4) y * 
/' \-'' *\ *V'Q'\,-'V\ 

/k 3 AX k s \-^-» g- £G(K + q) = G(K) K Q' K (.» 

* ** * * ■ 4 * » * (. 5i » >, 

K-Q'-Q 

FIG. 4. Use of the DCA Laue function Adca leads to a 
leads to the replacement of the lattice propagators G(ki), 
G(k.2), ... by coarse grained propagators G(K), G(K'), ... 

If we employ the DCA Laue function in each of the 
self energy diagrams shown in Fig. [I] then we may freely 
sum over the momenta within each coarse-graining cell 
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shown in Fig. 0. As illustrated in Fig. [| for a fourth- 
order graph, this leads to the replacement of the lattice 
propagators G(ki), G(k 2 ), ... by coarse grained prop- 
agators G(K), G(K'), .... which are given by 



(10) 



where N is the number of points of the lattice, N c is 
the number of cells in the cluster, and the k summation 
runs over the momenta of the cell about the cluster 
momentum K (cf. Fig. |^). 

As N c increases from one, systematic non-local cor- 
rections to the CPA self energy are introduced. To see 
this, recall that the self energy is a functional of the 
Green function. According to Nyquist's sampling the- 
orem to reproduce correlations of length < L/2 in 
the Green function and corresponding self energy, we 
only need to sample the reciprocal space at intervals of 
Afc w 2ir/L. Knowledge of these Green functions on a 
finer scale in momentum is than AA; unnecessary, and 
may be discarded to reduce the complexity of the prob- 
lem. Thus the cluster self energy will be constructed 
from the coarse-grained average of the single-particle 
Green function within the cell centered on the cluster 
momenta. For short distances r < L/2, where L is now 
the linear size of the cluster, the Fourier transform of 
the Green function G(r) « G(r) + 0((rAfc) 2 ), so that 
short ranged correlations are reflected in the irreducible 
quantities constructed from G; whereas, longer ranged 
correlations r > L/2 are cut off by the finite size of the 
cluster §. 

We show in the appendix that that free-energy ar- 
guments presented previously j9) apply to the disor- 
dered case as well. In particular, the DCA estimate 
of the lattice free energy is minimized by the choice 
E(k,w)=E(M(k),w). 

Algorithm With the substitution A — * Adca, 
most of the diagrams represented in Fig. |l| remain. 
However, the complexity of the problem is greatly re- 
duced since the nontrivial sums involve only the cluster 
momenta K ( numbering jV c instead of N). Further- 
more, since these diagrams are the same as those from 
a finite-sized periodic cluster of N c sites, we can eas- 
ily sum this series to all orders by numerically solving 
the corresponding cluster problem. The resulting al- 
gorithm is as follows (for notational convenience the 
frequency arguments are not displayed below): 

1. Make a guess for the impurity- averaged cluster 
self energy E(K), usually zero. 

2. Calculate the coarse-grained Cluster Green func- 
tions 



«4e 



N i u + M ~ " K + k ~ S(K) 



(ii) 



3. Calculate the cluster-excluded 
propagator Q{K) = 1/(1/G(K) + E(K)). The 
introduction of G(K.) is necessary to avoid over- 
counting diagrams on the cluster. 

4. Fourier transform Q to the real space matrix rep- 
resentation of the cluster problem, (i.e., write 
Gn,m = Hk <?(K)expiK • (r„ - r m ) ) whence 
the disorder potential may be represented as a 
diagonal matrix V (with elements V n on the clus- 
ter sites labelled by n) and form a new estimate 
of the disorder-averaged cluster Green function 
matrix 



(12) 



where the average () indicates an average over 
disorder configurations on the cluster Q . 

5. Transform back to the cluster reciprocal space 
and form a new estimate of the self energy 
E(K) = 1/G(K) - 1/G(K) 

6. Repeat, starting from 2, until E(K) converges to 
the desired accuracy. 

The algorithm recovers the CPA for N c = 1 and be- 
comes exact when N r — > oo. 



IV. CHARACTERISTICS OF THE DCA 

In Ref. H A. Gonis discusses the CPA, and vari- 
ous methods to incorporate non-local corrections. He 
lists the most important characteristics of a success- 
ful cluster theory. A successful theory must be able to 
account for fluctuations in the local environment in a 
self-consistent way, become exact in the limit of large 
cluster sizes, and recover the CPA when the cluster size 
becomes one. It must be straight-forward to implement 
numerically and preserve the full point-group symme- 
try of the lattice. Finally, and most significantly, it 
should be fully causal so that the single-particle Green 
function and self energy are analytic in the upper half 
plane. In the next two sections, we demonstrate that 
the DCA satisfies each of these requirements. 

The limits N c — > 1 and N c — > oo : As mentioned 
above, the DCA recovers the CPA for N c = 1. When 
N c = 1, K = 0, and k = k. Then the DCA algorithm 
reduces to the self consistent scalar equations (in con- 
trast to the DCA which involves matrix equations): 



G 



N 4-^ 



1 



LU + fl - £k 



E ' 



g- 1 = 1/g + e, 



(13) 



(14) 
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G=((g- 1 -V)- 1 ), (15) 

which together correspond exactly to the prescription 
for CPA H , and are easily solved, for example by iter- 
ation. 

As N c — > oo the DCA becomes exact, including 
correlations over all length scales. For, the DCA 
Laue function requires complete momentum conserva- 
tion in this limit, i.e., [K,G(K)] = [k, G(k)], whence 
Q = l/(cj + /i — £k) which is the bare lattice Green func- 
tion so that Eq. 12 amounts to solving the problem by 
exact diagonalization. 

Causality : It is easy to show that this algorithm 
is fully causal || . 
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FIG. 5. Illustration of the essential steps of the proof 
that the DCA is causal (see text). 

For our purposes, this is equivalent to requiring that 
all the retarded propagators and self energies are Her- 
glotz, or analytic in the upper half plane. Since the 
diagrams which describe the DCA cluster problem are 
isomorphic to those of a real finite-size periodic clus- 
ter, the corresponding impurity averaged cluster self 
energy shares the causal properties of this system. Fur- 
thermore, the coarse graining step (2) cannot violate 
causality, since the sum of analytic functions is ana- 
lytic. The only "suspect" step is the cluster-exclusion 
step (3), thus we must show that the imaginary part 
of Q is negative semidefinite. The essential steps of the 
argument are sketched in Fig. |^. The imaginary part of 
G(K,cu) = (GiK^uj)- 1 + E(K,u))-} is negative pro- 
vided that Im^K,^)- 1 ) > -Im£(K,w). G(K» 
can be written as G(K, uj) — (N c /N) X)k( z K+k) _1 ( aj )' 
where the z K+k (w) are complex numbers with a posi- 
tive semidefinite imaginary part — Iml](K,cj). For any 
K and uj, the set of points z K+k (w) are on a segment of 
the dashed horizontal line in the upper half plane due 
to the fact that the imaginary part is independent of 
k. The mapping z — > 1/z maps this line segment onto 
a segment of the dashed circle shown in the lower half 
plane. G(K, uj) is obtained by summing the points on 



the circle segment, yielding the empty dot that must lie 
within the dashed circle. The inverse necessary to take 
G(K,w) to l/G(K,u;) maps this point onto the empty 
dot in the upper half plane which must lie above the 
dashed line. Thus, the imaginary part of G(K,o;) _1 is 
greater than or equal to — ImE(K,w). This argument 
may easily be extended for G(z) for any z in the upper 
half plane. Thus Q is completely analytic in the upper 
half plane. 

Preservation of the Lattice Symmetry : Since the 
DCA is formulated in reciprocal space, it preserves the 
translational symmetry of the system. However, care 
must be taken when selecting the coarse-graining cells 
to preserve the point group symmetries of the lattice. 
For example, the calculations presented in Sec. [v| are 
done for a simple square lattice. Both it and its recipro- 
cal lattice have a C<± v symmetry with eight point group 
operations. We must choose a set of coarse-graining 
cells which preserve this point-group symmetry. This 
may be done by tiling the real lattice with squares, and 
using the K points that correspond to the reciprocal 
space of the tiling centers. For large V, an important 
configuration of the half filled system is that with all 
of the sites with Vi — — V occupied and those with 
Vi = +V unoccupied. To retain this configuration on 
the cluster, when iV c > 1 we will choose N c even. 



o 


o 


o 


o 




/ o\ 


o 


o o 


o 


o 


o 


o 


o 


o 




/ o 


o 


cr 


\p/a 




o 


o 




Nc 


=4 




Nc=8 












o 


o 


o 


cN 


xO 


o 


°/ 


/o o 




X) 


o 


o 












o 


o o 


o 


(T 


















Nc= 


18 








o 


o 


o \ 


o 


o\ 


o 


o o 


o 


°/ 


/o 




Nc=10 
















cr 


\° 


o 


o 




o 


cA 


\p o 


°/ 


O 


o 


o 








o 


o 


o 




Xo 


O 


o 


o 


o 


o 


o 


o 




o 


o o 


\ 


O 


o 


o 


o 


o 


o 


o 




o 


o o 


o 


V) 


o 






Nc= 


=16 








Nc=20 








o 


o 


o 


o 


o 


(\ 


o 


o o 


o 


o\ 


o 


o 


o 


o 


o 


o 


o 




o o 






o 












a 
























\ 










o 


o 


o 


o 


o 


o 


cX 




o 


o 


o 



FIG. 6. Different tile sizes and orientations. The tiling 
principal translation vectors, ai and a2, form two sides of 
each tiling square (illustrated for the iV c = 20 tiling). For 
square tile geometries, a2 X = — a-iy and a2 y = a± x . 



Square tilings with an even number of sites include 
N c = 4, 8, 10, 16, 18, 20, 26, 32, 34, 36, • • ■. The first few 
are illustrated in Fig. ^. The relation between the 
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principal lattice vectors of the lattice centers, ai and 
&2, and the reciprocal lattice takes the usual form 
gi = 27raj/|ai x a 2 |, with K mn = ngi + mg 2 for 
integer n and m. For tilings with either a\ x = a\ y 
(corresponding to N c — 8, 18,32---) or one of a\ x 
or a\ y zero (corresponding to N c = 1, 4, 16, 36 • ■ •), 
the principal reciprocal lattice vectors of the coarse- 
grained system either point along the same directions 
as the principal reciprocal lattice vectors of the real 
system or are rotated from them by 7r/4. As a result 
equivalent momenta k are always mapped to equivalent 
coarse-grained momenta K. An example for N c = 8 is 
shown in Fig. @. However, for N c = 10, 20, 26, 34 • • •, 
the principal reciprocal lattice vectors of the coarse- 
grained system do not point along a high symmetry 
direction of the real lattice. Since all points within a 
coarse-graining cell are mapped to its center K, this 
means that these coarse graining choices violate the 
point group symmetry of the real system. This is illus- 
trate for N c = 10 in Fig. 0, where the two open dots, 
resting at equivalent points in the real lattice, fall in 
inequivalcnt coarse graining cells and so are mapped to 
inequivalcnt K points. Thus the tilings corresponding 
to N c = 10, 20, 26, 34 • • • violate the point-group sym- 
metry of the real lattice system and should be avoided. 



fective Boltzmann factor determined by the electronic 
partition function for that particular disorder configu- 
ration ||. 

In this section, we propose an approach to carrying 
out the disorder averaging by statistically sampling dis- 
order configurations using a Markov process. For sys- 
tems with quenched disorder, one could also sample 
random configurations of the disorder potentials and 
calculate the corresponding Green function using ma- 
trix inversion. A significant advantage of the Markov 
technique is that it may be easily modified to treat 
either quenched or annealed disorder |9). 

In a Markov process, each disorder configuration de- 
pends upon the previous configurations. To evolve 
from one configuration to another, we will propose lo- 
cal changes in the disorder potentials. These changes 
are accepted with some probability determined by ei- 
ther the effective Boltzmann factor ||, for annealed 
disorder, or the probability distribution P(Vi) , for 
quenched disorder. If we accept such a change in one 
of the disorder potentials, say on site I, then the new 
Green function matrix G' n rn depends on the previous 
Green function matrix G„ im through the matrix rela- 
tionship (where, once again, for notational convenience 
the frequency arguments are not displayed) 
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FIG. 7. The coarse graining cells for N c — 8 and 10 each 
centered on a coarse-grained momenta K represented as 
black filled dots. For N c — 8 equivalent momenta k are 
always mapped to equivalent coarse-grained momenta K. 
However, this is not true for N c — 10 where, for exam- 
ple, the two equivalent momenta shown by open dots are 
mapped to inequivalent coarse-grained momenta. 

An efficient numerical algorithm for disorder averag- 
ing : The implementation of the DCA clearly requires 
an efficient algorithm for disorder averaging on a clus- 
ter of size N c .( Needless to say, this particular aspect is 
common to all approaches where disorder averaging is 
involved.) Even for a system with binary diagonal dis- 
order, where each site can acquire only one of two val- 
ues for the potential (±V) the total number of disorder 
configurations is 2 °, which grows exponentially with 
N c . For a generic quenched disorder, the probability 
of the various disorder configurations are determined 
by the specified P(V); whereas, for annealed disorder, 
the probability of a configuration depends upon an ef- 
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Since we change on the potential on the site I and the 
matrices V and V are diagonal, their difference 5V = 
V — V is a diagonal matrix with only the Z'th diagonal 
clement finite. Then 
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If we set I = n in Eq. 17, we get 

Gi m 



Gj 



1 - Gi,i5Vi 

If we substitute this result back into Eq. |l7], we get 
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(19) 



an equation which requires roughly N% operations to 
evaluate |L4j for each frequency. 

Because the technique involves importance sam- 
pling, it is likely to miss rare configurations of disorder, 
and any special physics that arises from such configu- 
rations, such as Lifshitz tails in the DOS. However, 
when these configurations are well known, we can eas- 
ily adapt this method to include them. This may be 
done by excluding them from the sampling, and then 
including the corresponding configurations in the sam- 
ple with the appropriate reweighting. 
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V. RESULTS 

A. Single-Particle Properties 

To illustrate the algorithm discussed above, we 
present calculations on a two-dimensional square lat- 
tice system with t = 0.25 and binary random disorder 
so that Vi — ±V with equal probability. No tricks are 
used to force the algorithm to remain causal such as 
renormalizing the spectrum or cutting off any negative 
tails of the density of states. 
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FIG. 8. The self energy N c = 1 and for N c = 16 for 
different values of K for V = 0.1 (left) and V = 0.5 (right). 
The thin lines show E(K,w) for K = (0,0) K = (0,tt), 
K = (71-/2,71-/2) and K = (tt/2, 0), plotted versus w for 
N c = 16. For V — 0.1 the self energy at the different 
values of K is essentially the same as the N c = 1 result. 
For V = 0.5 they differ significantly, indicating that the 
momentum dependence of the self energy increases with V. 

The self energy is plotted in Fig. |8| for N c = 1 and 
for N c = 16 for different values of K for V = 0.1 (left) 
and V = 0.5 (right). For V — 0.1 the self energy for 
N c = 16 has very little momentum dependence, thus 
the different curves fall atop of one another. They also 
are hence very close to the self energy for N c — 1 (the 
CPA result) indicating that it is a very good approx- 
imation for the self energy when V is small. How- 
ever, for larger V = 0.5 the self energy curves at dif- 
ferent values of K for N c — 16 differ considerably 
from each other and from the CPA self energy obtained 
with N c = 1. Thus, as V increases non-local correc- 
tions clearly become important and are expressed in 
the momentum-dependence of the self energy. 
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FIG. 9. The density of states for N c = l and for iV c = 32 
for four values of disorder potential V. 

As shown in Fig. || (top left) the single particle 
density of states is essentially independent of N c for 
V = 0.1; however, for V = 1.0, the density of states de- 
pends strongly upon N c . In particular the gap around 
u> = 0, which is sharp for N c = 1, is partially filled in. 
The top and bottom of the band also acquire tails as 
N c increases. When N c > 1, the density of states ac- 
quires several additional structures which correspond 
to important local configurations of the disorder. The 
additional features and the band tails are absent in 
the CPA and believed to be due to local order in the 
environment of each site H . 



B. (Absence of) Localization 

Despite its advantages over the CPA as discussed 
above, one feature that the DCA shares with the CPA 
and similar self-consistent cluster methods is its limited 
ability to take into account localization effects |l^ [TtJ . 
To show this, we measure the probability that an elec- 
tron remains at site I for all time flla,[16fl: 



P(oo)= lim (\G(l,l,tf) 



lim^ / de (\G u {e + tv)' 2 



(20) 



As shown in |T^,[l(|], P(oo) is expected to be nonzero 
as long as there are a thermodynamically significant 
fraction of localized states in the spectrum of eigen- 
states of the disordered system. In 1 and 2 dimensions 
this is expected to happen for arbitrarily small but 
thermodynamically significant disorder. Since the clus- 
ter is formed by coarse-graining the real-lattice prob- 
lem in reciprocal space, local quantities on the cluster 
and the real lattice correspond one-to-one. Thus, to 
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test for localization, we need only apply the formula 
|20| for each site on the cluster. Making this substi- 
tution and introducing the local coarse-grained (but 
not disorder averaged) spectral function, A(l, uS) = 
— ilmCrj Eq. |2fl becomes 



P(oo) 



= lim p(r)) 

= lim 

rj-o N c 



FIG. 10. The local hybridization rate when V = 0.4 for 
several values of N c . The probability that the electron re- 
mains localized p(r\ — > 0), when e = (half filling) is shown 
in the inset. As r\ — > 0, p(r/) extrapolates to zero indicating 
the lack of localization. 

As described in Ref . , the hybridization rate be- 
tween the cluster and its host is given by 



E 



dtuduj' 



A(l,u)A(l,u') 
lu — lu' — 2irj 



(21) 



Im 



pirj) is plotted versus r\ in the inset to Fig. |l0| for the 
half-filled model when V = 0.4. The p(rf) extrapolates 
to zero, indicating the lack of localization. 

This result can be understood from either a diagram- 
matic perspective or by carefully assessing the cluster 
problem. As is well known []l7f , the crossing diagrams, 
especially those which involve many crossings, describe 
the coherent backscattering of electrons which are re- 
sponsible for localization. Hence the CPA, which in- 
cludes only non-crossing diagrams, can not describe 
localization. Within the DCA, however, for N c > 1 
some crossing graphs are restored . Within each dia- 
gram, each X represents scattering from distinct site. 
Since there are only N c sites on the cluster, the maxi- 
mally crossed DCA graphs can have at most N c cross- 
ings. Since all states are expected to be localized in two 
dimensional disordered system, apparently an infinite 
number of crossings are needed to describe localization 
diagrammatically. From the perspective of the clus- 
ter, this result is not surprising since each site on the 
cluster is coupled to a non-interacting translationally 
invariant host into which electrons can escape. Thus, if 
the density of states is finite at some energy, then the 
corresponding states can not be localized unless the 
hybridization rate at that energy between the cluster 
and the host vanishes. 



0.4 



0.3 



V=0.4 



3e-06 



2e-06 - 




§, 0.2 



\ le-06 - 

\ o r ■ < ■ 

\ 0.01 0.02 

\ 1 



\ 



r((»)N=i 

1 r(<D)N=8 

r(co)N=16 

r(o>)N C =32 




\G(K,l 



e(k, w ; 



(22) 



The net hybridization rate to a site on the cluster (the 
K-integrated r(K, lj)) is plotted in Fig. |l^ when V — 
0.4 for several values of N c . It remains finite over the 
entire region where the corresponding density of states, 
shown in Fig. is finite. This is consistent with the 
lack of localization demonstrated in the inset. 

We note that the hybridization falls as A^ c increases 
(for large N c , T(K, lj) ~ 0(1/N C ) @) especially at the 
band edges, although, given that it is defined entirely 
in terms of disorder-averaged propagators, it is still 
unlikely to be sensitive to localization effects. How- 
ever, the number of diagrammatic crossings in two par- 
ticle properties (such as the conductivity) which are 
strongly affected by localization effects does increase 
with N c even within the DCA. Thus, it is likely that 
disordered DCA can describe the precursor effects of 
localization. Some evidence for this can be seen in the 
(finite time) probability that an electron on a site I re- 
mains after a time t, P(t) = (\G(l, I, t)| 2 ^. As shown 

in Fig. [ll], for N c = 1, this probability falls quickly 
with time. The long time behavior is shown in the in- 
set. As N c increases, the electron remains localized for 
longer times. Hence one can hope that a careful finite 
size scaling study of two particle properties within the 
disordered DCA can even capture some aspects of the 
localization transition. 
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FIG. 11. The probability that an electron on a site I re- 
mains after a time t, P(t) = (\G(l, I, i)| 2 ) for several values 
of iV c when V = 0.4. 



VI. CONCLUSION 

We have developed a modification of the Dynami- 
cal Cluster Approximation to treat disordered systems. 
This formalism satisfies all of the characteristics of 
a successful cluster approximation. It is causal, pre- 
serves the point-group and translational symmetry of 
the original lattice, recovers the CPA when the cluster 
size goes to one, and becomes exact as N c — > oo. Like 
the CPA the problem is disorder averaged and has a 
simple diagrammatic formulation. It is easy to imple- 
ment numerically and restores sharp features and band 
tailing in the DOS which reflect correlations in the lo- 
cal environment of each site. Although the DCA does 
not capture the localization transition, it does describe 
the precursor effects of localization. It systematically 
restore the crossing graphs known to be responsible for 
localization, and might be able to access the localiza- 
tion transition itself via an appropriate finite size scal- 
ing analysis of two particle properties which remains 
to be developed. 

The DCA formalism we have discussed here can also 
be extended to problems with disorder and interactions 
simply by incorporating interaction diagrams in the self 
energy. This is also discussed in the Appendix below. 
The DCA should be able to provide a good descrip- 
tion of localization effects at finite temperatures in such 
contexts. For, in such cases the scattering processes are 
partially inelastic, so that the coherent back scatter- 
ing disappears after a characteristic inelastic scattering 
time |lq] . In this time only a finite- number of back- 
scattering processes can occur so only a finite number 
of diagrammatic crossings are needed to describe the 
finite-temperature physics, and these are captured in 
the DCA. 
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APPENDIX A: DISORDER-DC A FROM THE 
REPLICA METHOD 

An alternate way of justifying the DCA in the con- 
text of disordered systems is to use the replica (or other 
such) trick for disorder averaging p9| . For, this maps 
the disorder averaged problem into what looks like an 
interacting problem, whence the DCA formalism de- 
veloped by us earlier j^,|9| can simply be transcribed 
for this case, to arrive at the appropriate self consis- 
tent cluster problem. For the effective cluster problem, 
the replica trick can be " un-done" , and we recover the 
algorithm presented earlier in this paper. The same 
procedure also works for problems involving both dis- 
order and interactions. We detail this below. 

As is well known, for problems involving quenched 
disorder, as for example corresponding to the Hamil- 
tonian: 



H — H + Hdis 



k.cr 



Hdis = Vji 



(Al) 
(A2) 

(A3) 



where £ k = £k — and Vi is the random potential dis- 
tributed according to a given probability distribution 
P(V), complications arise because the disorder averag- 
ing has to be done on the free energy 

F = -k B T\nZ 

and the green functions 

Gij(r) = -Tr[T T C^(T)C}. a ex P (-/3H)]/Z. 

Here Z = Tr[exp(— f3H)\ is the partition function, and 
T T represents the imaginary-time ordering operator. 
In the replica trick (l9| , one writes 



In Z 



lim 



1 



and 1/Z = lim Z r 



and assumes that the order of taking the limit m r — > 
and disorder averaging can be interchanged. Then for 
any positive integer m r , the resulting disorder averaged 
quantities such as (Z mr ) , (G k (T)), etc., can be repre- 
sented in terms of an interacting problem involving m r 
replicas of the original electronic degrees of freedom, 
which we index with the subscript a = 1, ...,m r . 

For example, using the standard Fcrmionic (Grass- 
mann variable) functional integrals [^0| to represent 
the traces above, we can write 

(Z m "} = J Dc*Dc exp [-09] (A4) 

(G k (r)) =- J DcTDc c k)<Tll (r) c k>CTil (0) exp[-/W>] (A5) 

Here 9 is an effective free energy functional which 
arises from the disorder averaging, and can be writ- 
ten as 



0¥ = V / dr ci^ a (T){d r + a)c k , CT .a(r) + V W{n % ) 

(A6) 



k,cr,a ' 



where, 



r/3 

/ n-i,o-.a(V) dr 



and 



exp[-W(ni)] = (exp (-VSfii)) = / dV.P^ex^-V^) 

In terms of the cumulants (V l ) c of the disorder distri- 
bution P(V), one can write 



1=2 
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So, clearly, W introduces (local in space but non-local 
in time) interactions between electrons belonging to 
arbitrary replicas. 

If one re-expands exp [— W(fij)] in powers of the cu- 
mulants (V*) one can perform the Fermionic traces 
using the standard techniques of diagrammatic pertur- 
bation theory. Then, order-by-order in perturbation 
theory, the dependence on m r is explicit and analytic, 
and the lim mr ._>o can be evaluated precisely. The re- 
sulting terms are in exact, one-to-one correspondence 
with the terms obtainable by writing out the diagrams 
from a direct perturbation expansion in powers of Vi 
and then disorder averaging as discussed in Section II. 
The m r — + limit eliminates the diagrams (in the in- 
teracting problem ) containing internal loops with free 
sums over the replica indices (as required, since such 
diagrams never appear in the direct disorder averaged 
perturbation theory formalism of Section II ). 

For the "replicated interacting problem" obtained 
above, one can transcribe exactly the DCA formalism 
discussed in refs. ]|||. If one assumes that the self- 
consistent host propagators do not break replica sym- 
metry, then the effective cluster problem corresponds 
to a Fermionic functional integral involving an effec- 
tive, self consistent cluster free energy functional given 

by 

r/3 rP 



tive free energy functionals for the lattice and the clus- 
ter pick up the additional terms U ^ Q J2i N n i,ha n i,l,a 
and U J2 a J2i N c n i,ha n i,l,a- The resulting cluster 
problem now has both interactions and disorder on 
the cluster of N c sites which are dual to the cluster 
momenta K: a bare ( retarded ) cluster propagator 
(/ _1 (K, t — t'), a random potential Vi distributed ac- 
cording to P(Vi), and the Hubbard interaction U at ev- 
ery site i of the cluster. One can resort to any technique 
of one's choice to solve this problem for the disorder- 
averaged cluster Green functions G(K, uj) and cluster 
self-energies S(K, w) and go through with the rest of 
the DCA iteration. 



^c=J2 I dT i dr 1 c^ a . a (T)Q-\K,T-T')cK,„. a {T') 

K.,a,at J ° 

+ Y,W(n t ) (A7) 



But, as is easy to see using the same procedure as out- 
lined earlier in this appendix, such an effective free en- 
ergy functional is exactly what one would obtain if one 
were to disorder average (using the replica trick) a clus- 
ter problem with N c sites which are dual to the cluster 
momenta K, a bare ( retarded ) cluster propagator 
(? _1 (K, t — t'), and a random potential Vi distributed 
according to P(VJ) at every site i of the cluster. Hence 
we have an alternate justification for the disorder-DCA 
algorithm set down in Section III. The above route also 
enables one to quickly extend our discussions in ref. || 
regarding the 2-particle propagators, Ward identities, 
etc., to the disorder-DCA context. Most significantly, 
the DCA estimate of the lattice self energy is mini- 
mized by the choice E a (k,w) = E(M(k), u>). 

We note that the arguments presented in the main 
text and in this appendix are also easily extended to 
problems involving interactions and disorder. For ex- 
ample, for the case of the Hubbard model with diagonal 
disorder, one would add to the starting Hamiltonian 
the interaction term U J2i n i,1 n i,l- Going through 
exactly the same procedures as outlined above, it is 
not hard to see that the only change is that the effec- 



11 



